Theory of Interferometric Photon-Correlation Measurements: Differentiating 

Coherent from Chaotic Light 
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Interferometric photon-correlation measurements, which correspond to the second-order intensity 
cross-correlations between the two output ports of an unbalanced Michelson interferometer, are 
sensitive to both amplitude and phase fluctuations of an incoming beam of light. Here, we present the 
theoretical framework behind these measurements and show that they can be used to unambiguously 
differentiate a coherent wave undergoing dynamical amplitude and phase fluctuations from a chaotic 
state of light. This technique may thus be used to characterize the output of nanolasers and monitor 
the onset of coherent emission. 
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I. INTRODUCTION 

The laser threshold highlights a radical change in the 
mode of operation of the device. When the laser operates 
at powers below threshold, the emission of the excited 
dipoles constituting the active medium is understood to 
be spontaneous (such that on average there is less than 
one photon in the cavity) and the laser output is inco- 
herent. The incoherent field is described as a statistical 
ensemble of waves emitted at random times with random 
phases and displays chaotic (thermal) photon statistics 
|l|, [1]. On the other hand, at powers above threshold, 
the emission is stimulated (as the number of photons in 
the cavity exceeds one) and the laser output is coher- 
ent. The state of the field is described by a quantum 
mechanical coherent state, which is usually represented 
in phase-space as a rotating vector (albeit with quantum 
uncertainty) called a "phasor" , and displays Poisson pho- 
ton statistics [I]- 

The coherence of the laser output is usually evaluated 
by measuring two types of autocorrelation functions. The 
first-order (field) autocorrelation function is usually mea- 
sured by means of an interferometer that introduces a 
delay r between its two arms, and is defined as 



{E*{t + r)E{t))t 
{E*{t)E{t))t 



(1) 



where E*{t)/ E{t) is the positive/negative frequency part 
of the electric field operator at time and {■.■)t denotes 
a quantum mechanical and statistical average, as well as 
an average over time t. For an ergodic stationary field, 
the statistical averaging is equivalent to the time-average, 
so that in that case the autocorrelation function may be 
evaluated at an arbitrary time, e.g. at t = 0. The first- 
order autocorrelation function is essentially the Fourier 
transform of the field spectrum, and thus for a Lorentzian 
spectrum it corresponds to a decaying exponential, 

9''Hr) 



— iujT — |r I /t 



(2) 



with Lorentzian spectra give first-order correlation func- 
tions of the form of Eq. (2). Measuremant of Tc permits 
nevetheless a quantitative distinction between the chaotic 
or coherent output of a laser, as the value of Tc becomes 
longer above threshold, reflecting the narrowing of the 
Schawlow-Townes linewidth of the laser as a function of 
pump power. However, as Glauber pointed out in his 
seminal papers on the theory of coherence: "Coherence 
does not require monochromaticity. Coherent fields can 
be generated with arbitrary spectra" implying that 
a measurement of g^^^ is insufficient for identifying a co- 
herent state. 

The second-order (intensity) autocorrelation function, 
defined as 



{E*{t)E*{t + T)E{t + T)E{t))t 



(3) 



where uj is the frequency of the field and Tc is its co- 
herence time. Both chaotic and coherent states of light 



is usually measured in a Hanbury Brown and Twiss 
photon-counting experiment, as the probability of de- 
tecting two photons separated by a time-difference r. 
The second-order autocorrelation function measures es- 
sentially the time-distribution of the relative intensity 
fluctuations and its value at zero time-difference (r = 0), 
gives the variance (cr^) of the relative intensity fluctua- 
tions of the field through the expression = f/^^^ (0) — 1. 
The value of g^^-'(O) is often used to characterize the co- 
herence of the laser output, as it is equal to 2 for chaotic 
light and equal to 1 for a stable coherent wave with no 
amplitude fiuctuations. 

In conventional (large) lasers, when the pump power is 
increased across the threshold, there is an abrupt transi- 
tion in the value of the zero time-difference second-order 
intensity autocorrelation from g^^\0) = 2 to g^^^{0) = 1, 
a feature that has led to the identification of this tran- 
sition with the threshold [s^. However, in nanolasers the 
transition between the two values of g^'^\0) as a func- 
tion of input power is very gradual, even when an abrupt 
change in output power clearly identifies a threshold. 
This feature is often interpreted as indicating that the 
output of nanolasers is not fully coherent above thresh- 
old, but includes incoherent spontaneous emission P-ffil- 
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This interpretation rests on the assumption that the 
value of g*-^^(0) has a direct (bijective) relationship with 
coherence. However, while the value ^'■^^(O) = 1 can 
only be obtained for a coherent stable wave, the val- 
ues g^^HO) > 1 do not necessarily correspond to the 
presence of incoherent chaotic light in the laser out- 
put. They may also arise from a coherent quantum 
state undergoing dynamical amplitude fluctuations. In 
this context, we note that in conventional large lasers 
the value ^'^^■'(O) = 1 arises because their gain (i.e. the 
number of excited dipoles) is clamped above threshold, 
so that the output field presents no amplitude fluctua- 
tions if the pumping rate is stable. In a laser, the co- 
herent field undergoes phase fiuctuations (corresponding 
to the Schawlow-Townes phase diffusion) and these give 
rise to its finite spectral width, while preserving the value 
g^'^\0) = 1 since the gain remains clamped. On the other 
hand, if the gain were not clamped but fluctuated in time, 
it would produce amplitude (intensity) fluctuations in the 
laser output, giving rise to g^^^(O) > 1. Thus, the value of 
g^^\0) alone is not sufficient for distinguishing between 
chaotic and coherent light. 

In a recent publication, we presented a novel exper- 
imental technique, consisting of interferometric photon 
correlation measurements, which gives qualitatively dif- 
ferent results when applied to chaotic and to coherent 
light and can thus unambiguously discriminate between 
these two types of states of the electromagnetic field 
This technique, a variant of which was initially devel- 
oped to study the dynamical fluctuations in the fluores- 
cence spectrum of single molecules [1, can measure 
the correlation times of both phase and amplitude fluc- 
tuations of an incoming beam. It can thus distinguish a 
coherent field from a statistical distribution of randomly- 
phased waves (chaotic field) as in the latter case the two 
types of fiuctuations have a fixed relationship. By use of 
this technique, we showed that the output of a nanoscale 
laser is indeed a coherent state undergoing amplitude 
fluctuations, and not a chaotic statistical ensemble of 
spontaneously-emitted waves. We thus conflrmed ear- 
lier publications in which we presented evidence that the 
fluctuations in the output of a nanolaser is not due to the 
presence of incoherent spontaneous emission in the out- 
put of the laser, but is the result of the laser's relaxation 
oscillations excited by the "discretization noise" associ- 
ated with the small (and integer) number of photons at 
threshold [10]. 

In this paper, we develop the theory of interferometric 
photon correlation experiments and examine the corre- 
sponding correlation function for three types of states 
of the electromagnetic held: (1) chaotic light (2) sta- 
ble coherent light and (3) coherent light with amplitude 
fluctuations. We show that the interferometric cross- 
correlation function for chaotic light is qualitatively dif- 
ferent from that of coherent light, thus identifying a 
signature for coherent light. We thus show that even 
though chaotic light and amplitude-modulated coherent 
light may have similar first- and second-order correlation 



functions, they may be distinguished by their interfero- 
metric cross-correlation functions. 



II. INTERFEROMETRIC PHOTON 
CORRELATION MEASUREMENTS 

We consider a light beam of central frequency w and 
coherence time Tc (for simplicity we take the spectrum 
to be Lorentzian), entering an unbalanced Michelson in- 
terferometer (see Fig. [1} in which one arm is longer by 
a distance d, thus introducing a delay of d/c (where c is 
the speed of light) with respect to the other. The elec- 
tric field operators at the two output ports A/B of the 
interferometer may be written in terms of the input field 
operators as 

E-M - ^; °'""''(|-°''" 

E-M = E.5!(i±i?Zf)±i5!W ,4) 

v2 

where Eq and a^{t) are, respectively, the vacuum field 
and the creation operator for the input mode. An equiv- 
alent pair of equations holds for Ea/Eb and the anni- 
hilation operators. Using these relations, the intensity 
coming out at each port for a stationary input field, may 
be written as 

{EimAit)) = |(l-25R[g(i)(d/c)]) 
{EUt)EBit)) = |(l + 25R[g(i)(d/c)]) (5) 

where Ii — \Eo\'^{a^a) is the input intensity and di[x] is 
the real part of the complex number x. 

The quantity that is measured in these experiments is 
the second-order intensity cross-correlation function be- 
tween the two output ports of the interferometer, which 
may be written as: 

ff(2^)(r,d/c) = {EmE*BiT)EB{r)EAm (6) 

where t is the time-difference between the photon detec- 
tion events in ports A and B, while the normalization 
of g^^'^^ through a denominator equal to the product of 
the two mean output intensities is implicit. Expressing 
terms of the input field operators through Eqs. 
16 terms are obtained. Among these terms, those 
that involve an unequal number of creation and annihi- 
lation operators delayed by d/c (for example, the term 
{a^O)a^d/c + T)a{T)aiO))) oscillate like e'""^/" or e^'""^/" 
(or their complex conjugates) and thus average out to 
zero when the interferometer arm difference d is dithered 
over a distance of a few wavelengths. The six remaining 
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FIG. 1: Experimental setup for interferometric photon corre- 
lation measurements. Light enters an unbalanced Michelson 
interferometer. The pathlength difference between the two 
arms of the interferometer d is dithered over several wave- 
lengths to average out the interference fringe contrast. Two 
single-photon detectors at the two output ports of the inter- 
ferometer (A/B) measure the cross-correlation in the photon 
arrival times. 



The first two terms correspond each to g^'^^r), the 
"standard" second-order autocorrelation function of the 
input field, and arise from the two photons propagating 
like particles along the same interferometer arm, while 
the third and fourth terms correspond respectively to 
g^^^T + d/c) and g^^^r ~ d/c) (that is to the second- 
order autocorrelation function displaced by r = —d/c 
and T = d/c) and arise from the two photons propagating 
like particles along different interferometer arms. These 
four terms are sensitive only to the intensity fluctuations 
of the field. The last two terms, which have a negative 
sign, describe the fluctuations of the field when it under- 
goes interference by propagating through both arms of 
the interferometer. It is these terms, which will be called 
hereafter "interference terms" , that permits us to differ- 
entiate chaotic from coherent light, as they are sensitive 
to both amplitude and phase fluctuations. They corre- 
spond to the conditional probability that a constructive 
interference of the two paths with difference d results in 
the detection of a second photon a.i t — t, given that a 
constructive interference at t = led to the detection of 
a first photon. Eq. (O may thus be re-written as 



terms are: 



g^^^\r,d/c)^U 



a\d/c)a\d/c + T)a{d/c + T)a{d/c) 

a\{))a\d/c + T)a{d/c + r)a(O) 
a\d/c)a\T)a{T)a{d/c) 
a){d/c)a){T)a{d/c + T)a{Q) 
a)(S))a){d/c + T)a{T)a{d/c)) (7) 



7(2^) (r, d/c) = i {2.9(2) (r) + {d/c + t)+ 5(2) {d/c - t) - 2^[g(^'^ {d/c, r, d/c + r, 0)] } 



(8) 



r 



where g^^-* (^4, is, ^2, ^i) = (aU^4)'^U^3)'^(^2)a(ii)) is the 
four-time second-order correlation function of the input 
field that corresponds to the "interference terms" . 

We note that when the interferometer is balanced 
{d — 0), the interference terms reduce to the standard 
intensity autocorrelation function and thus the interfer- 
ometric cross-correlation function reduces to 



9^'''\r,0) = y^-'>{r) 



1 



(2)^ 



(9) 



where the factor 1 /2 arises from the averaging over sev- 
eral interference fringes due to the dithering: When the 
interferometer is exactly balanced, light comes out only 
at port B and the cross-correlation is equal to zero; how- 
ever, when the interferometer is between two fringes light 



comes out from both ports, corresponding to a cross- 
correlation of one, giving an average of 1/2 upon dither- 
ing. However, when the interferometer is unbalanced, in- 
troducing a delay longer than the coherence time of the 
incoming radiation d/c > Tc, the interferometric cross- 
correlation function deviates from Eq. (9) and its fea- 
tures depend on the nature of the incoming radiation. 
In the next few sections, we will examine the features of 
g^'^^^T,d/c) for a chaotic field, described as a statistical 
ensemble of waves and for a coherent field, described ei- 
ther as a single stable wave with phase fluctuations, or 
as a wave with both phase and amplitude fluctuations. 



4 



III. CHAOTIC FIELD 

We consider the electromagnetic field produced by a 
laser pumped below threshold, in which TV excited dipoles 
emit spontaneously and independenty, each one photon. 
This field may be described as a statistical ensemble of 
N randomly-phased waves, so that its electric field may 
be written as 



N 



(10) 



where the phase 4>j , associated with the wave emanating 
from emitter j, undergoes random jumps with correlation 
time Tc, according to the equation of motion. 



(11) 



where randomly fluctuating frequency fj is a Langevin 
noise source obeying the standard correlation relations, 

im) = 

(/,(i2)/fe(<i)) = -SjkS{t2^h) (12) 



where Sjk is the Kronecker (5-symbol, while S{t2 — ti) is 
the Dirac (5-function. 

The first-order correlation function of this field may be 
calculated as 



N 



N 



)-#fc(0) 



g-«"-rg-kl/Tc 



(13) 



while the second-order correlation function is: 



N 



J<pj{0)+i4>kiT)-i<pi{T)-i(t>^{Q)\ 



j.k.i.m 

= l + |g(i)(T)|2 = l + e-2kl/-c 



(14) 



It corresponds to a peak at r = with a value of 5'^^ (0) = 
2 and decaying to g^^^(oo) — 1 exponentially at a rate 
2/tc. The interference terms in g^^^^ (Eq. [8|) are: 



1 ^ 

I 



(15) 



SO that the interferometric second-order intensity cross- correlation function g^^'^' for a chaotic field is, 
I 



(16) 



r 



Thus, it presents two peaks at r = ±c?/c (see Fig. [2]), 
corresponding to a pair of initially bunched photons trav- 
eling each on a different arm of the interferometer. Its 
value at T = is g'''^\Q,d/c) = 1, independent of the 
path-length difference d (averaged over several interfer- 
ence fringes), since the two photons of two interfering 
elementary randomly-phased waves have equal probabil- 
ity of coming out in either one of the two output ports. 



IV. STABLE COHERENT FIELD 



magnetic field can be described as a coherent quantum 
mechanical state that consists of a superposition of all 
number-states of the lasing mode, such that it presents a 
non-zero and oscillating expectation value for the electric 
field operator. The expectation value of the electric field 
thus corresponds to a classical oscillating stable wave. At 
the same time, the laser field undergoes phase diffusion 
describable by the Schawlow-Townes process, so that the 
electric field may be written as 



E*{t) = i;*e*"*+*'^(*) 



(17) 



We consider the light produced by a conventional laser 
operating above threshold. The state of the electro- 



where = Eq^J (at a) is the coherent amplitude of the 
electric field, while the phase undergoes random jumps 
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FIG. 2: Interferometric second-order intensity cross- 
correlation function for a chaotic field gl^^otici''~y'^/^) 
output ports of an unbalanced Michelson interferometer with 
pathlength difference d = 5crc (continuous red curve) and 
d = lOcTc (dashed green curve). 

with correlation time given by the Schawlow-Townes 
phase diffusion time according to 

'i^ = U (18) 

where is a randomly fluctuating Langevin noise source 
obeying 

ifm = 

{Uit2)Uiti)) = -Sih~h) (19) 
The first-order correlation function of this field is 

= e-^"^e-l^l/^* (20) 

the second-order correlation function is 

g(2)(T) = 1 (21) 

reflecting the lack of photon bunching in a coherent state. 
The interference terms in g^^^^ (Eq.[S]) may be calculated 
as: 

g^^^id/c,T,d/c + T,0) = (e'0('i/c)+«0(r)-»0(d/c+r)-»0(O)^ 

= e"2"'/^* (22) 

where m — min{\d/c\, \t\) is the smaller of the two values 
\d/c\ or |t|. Thus the interferometric second-order inten- 
sity cross-correlation function for a stable coherent field 
undergoing phase diffusion is 

5^of (^,rf/c) = l-^e-2Wr. (23) 

It consists of a flat baseline at g'-'^^\oo,d/c) = 1 — 
ig-2d/(cT0)^ in which there is a dip at r = going down 



to the value of g^^^\0,d/c) = 0.5 (see Fig. [S]). This 
dip is the signature of a coherent state and is most pro- 
nounced when the interferometer is strongly unbalanced 
{d ^ CT^) as in this case the baseline assumes its highest 
value g^^-^\oo,d/c) = 1. Its full- width measured at the 
value of 1 — l/(2e) — 0.816 gives directly the correlation 
time of the random phase jumps, r^. Its value at r = 0, 
significantly below 1, indicates that if a first photon is 
detected in one particular output port (say, in port A), 
a second photon arriving at times t <^ will also exit 
by the same port since for those short times the coherent 
state retains the memory of its phase and constructive 
interference will occur for the same output port. At long 
times however, t 3> t^, the second photon may exit from 
either port. 
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FIG. 3: Interferometric second-order intensity cross- 
correlation function for a stable coherent wave g^^^\T,d/c) 
at the output ports of an unbalanced Michelson interferome- 
ter with pathlength difference d — O.Scr^ (dotted red curve), 
d — CTcj, (dashed green curve) and d ^ ctc (continuous black 
curve) . 



V. COHERENT FIELD WITH AMPLITUDE 
FLUCTUATIONS 

A coherent field undergoing strong amplitude fluctu- 
ations in addition to Schawlow-Townes phase diffusion 
may be described by a simple phcnomenological model 
consisting of a wave modulated in amplitude (or in inten- 
sity) by an external signal whose frequency fmod varies 
randomly, following the correlation relations 

ifmodit)) = 
{Uod{t2)fmod{tl)) - :-^S{t2-ti) (24) 

where Tint is the correlation time of the intensity fluctu- 
ations. We may further assume for simplicity that the 
amplitude and phase fluctuations are uncorrelated, and 
that the amplitude modulation depth is large, so that the 
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unmodulated part of the carrier wave may be ignored. 
The electric field may thus be written as, 



E*{t) ^ V2i;*cos(x(t))e*"*+*'^(*) 



(25) 



where the modulation phase x(t) is driven by the 
randomly- varying frequency according to 

m = fmod (26) 

The first-order correlation function is 
5(i'(t) - 2(cos(x(r))cos(x(0)))e--(e^^M-^'^(")) 



gi"i-g-|i-|/(4T„t)g- 



where the correlation time Tc includes both phase and 
amplitude (intensity) fluctuations, 



1 _ f ^ 1 



(28) 



The second-order correlation function for a wave un- 
dergoing amplitude fluctuations may be calculated as 



g('\T) = 4(cos2(x(r))cos2(x(0))) 
= f 



(29) 



It thus resembles the second-order correlation function 
for a chaotic field with a peak at r = 0, but unlike 
the chaotic field that peak is not related to the first- 
order correlation function g^^^^ (r) and the spectrum of 
the field. It is nevertheless related to its power spec- 
trum, so that Eq. (29) may also be expressed in terms 
of the Fourier transform of the Relative Intensity Noise 
(RIN) spectrum, FT[RIN], as 



c,(2)(^) ^ l + FT[RIN] 



(30) 



which, in turn, is defined in terms of the spectral density 
of the intensity fluctuations, Ssi, as 



RIN = 



11 



(31) 



The interference terms in (^(^x) [g]) u^ay be calculated 



5^2) (d/c, r, djc + r, 0) = 4(cos(x(d/c)) cos(x(t)) cos(x(d/c + r)) cos(x(0)))(e''^(''/^)+'^(")-''^('^/^+")-'^(")) 

1+ le-MM„,^g-2mA 



(32) 



r 



where M — max(|d/c|, |t|) is the larger of the two values jd/cj or |r|, so that the g'^'^'^^ may be written as, 
I 



where the terms that give rise to features in the central 
part of the curve (i.e. near r « 0) have been grouped 
inside square brackets [...] for clarity. 

The interferometric second-order intensity cross- 
correlation function for a coherent field undergoing am- 
plitude modulation at a randomly varying frequency (see 
Fig mil bears the signature of a coherent state which is a 
dip at r w going down to 5i?^-//„ct(0' '^/c) = 0.75 and 
displaying a width of 2/tc. However, in contrast to the 
case of a stable coherent wave, the presence of amplitude 
fluctuations cause the central part of 9amv-fiuct^'^^'ll'^) 



to rise to a broad peak of width 1/Tmt inside which the 
sharp coherence dip appears, since 2/tc < l/Tint- Thus, 

the central part of gamp-fiucti''''^/'^) gives both the co- 
herence time of the laser and the correlation time of its 
amplitude fluctuations. Away from the central region, 

(2X) 

9amp~fiuct^^-< '^/^) displays two peaks at r = ±d/c corre- 
sponding to two displaced replicas of the autocorrelation 
function, as is also the case for the chaotic field. 

It should be noted that a "mixture" consisting of a 
fraction x coherent and (1 — x) chaotic light displays an 
interferometric cross-correlation function of the form 
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,-2|<i/c+T|/re ^ 



1 



,-1\dlc-r\lr^ _ 



= -2M|/c 



(34) 



where the term that gives rise to a feature in the central 
part of the curve (i.e. near r w 0) is inside square brack- 
ets [...] for comparison with Eq. As can be seen on 
Fig- El for a coherent-chaotic "mixture" the coherence dip 
appears on a flat baseline whose value is 1, for d 3> cr^. 
There is no peak rising above the value of 1 in the cen- 
tral part of this curve, in contrast to the case for coherent 
light with amplitude fluctuations (compare with Fig.|4]). 
Thus, the presence of the broad central peak together 
with a dip in the interferometric cross-correlation func- 
tion (arising from the terms in square brackets in Eq. 
(j33p ) permits to unequivocally discriminate a coherent 
field undergoing amplitude fluctuations from a "mixture" 
of coherent and chaotic light, even though both types of 
fields display similar second-order correlation functions 




-150 -100 -50 50 100 150 
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FIG. 4: Interferometric second-order intensity cross- 
correlation function for a coherent field with intensity fluc- 
tuations gf'^p- fiucti'^i'^/ output ports of an unbal- 
anced Michelson interferometer with pathlength difference 
d = SOcTc = 5cTint (continuous red curve), d = IOOctc = 
lOcTint (dashed green curve). The characteristic time of the 
intensity fluctuations is taken to be Tint = Wtc- 



VI. SUMMARY AND CONCLUSION 

Interferometric photon correlation measurements, in 
which the beamsplitter of a standard Hanbury Brown and 
Twiss setup is replaced by an unbalanced Michelson in- 
terferometer, are sensitive to the dynamical fluctuations 
of the incoming field and thus can discriminate between 
a field fiuctuating in time and a statistical ensemble of 
randomly-phased waves. They can therefore be used to 
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FIG. 5: Interferometric second-order intensity cross- 
correlation function for a 50-50 "mixture" of coherent and 
chaotic light gl^i^\d/c,T) at the output ports of an unbal- 
anced Michelson interferometer with pathlength difference 
d = 5cTc (continuous red curve), d = IOctc (dashed green 
curve). 



differentiate coherent (but fluctuating) light from chaotic 
light and thus pinpoint the onset of coherence of a laser. 

The unbalanced Michelson interferometer, which su- 
perposes two time-separated parts of the fleld, embodies 
Glauber's pioneering view of coherence: "In physical op- 
tics the term is used to denote a tendency of two values of 
the field at distantly separated points or at greatly sepa- 
rated times to take on correlated values... The coherence 
conditions restrict randomness of the fields rather than 
their bandwidth" [lt|. Alternatively, the measurement of 
the second-order intensity cross-correlation function at 
the output ports of the interferometer is equivalent to a 
photon-counting self-homodyne detection in the time do- 
main, with the signal rcfiected back from the long arm 
corresponding to the local oscillator. 

The second-order cross-correlation function at the 
two outputs of an unbalanced Michelson interferometer, 
g^'^^^T, d/c), is qualitatively different for chaotic and for 
coherent fields. While gchaoud'^^ '^/'^) ^ chaotic field 
consists of two peaks corresponding to two replicas of g'^^ 
displaced to r = ±d/c and is fiat around t w assum- 
ing the value g'^^^^O, d/c) = 1, the second-order cross- 
correlation function for a coherent field with amplitude 
fluctuations, in addition to the two replicas of g^^-' dis- 
placed to r = ±d/c, it displays a dip at t = whose 
width is given by twice the inverse of the coherence time 
of the field, superimposed on a broad peak that corre- 
sponds to g^^^r). These features, consisting of a sharp 
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dip dug into a broad peak, constitute the signature of 
a coherent field undergoing amphtude fluctuations and 
are absent for a chaotic fleld. Thus, inspection of the 
central part of g^^-^^ (r, d/c) (i.e. around r « 0) permits 
an unequivocal discrimination between a coherent and a 
partially chaotic field. 
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